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ABSTRACT 

We present a new calibration method based on cross-correlations with WMAP and apply it to 
data from the Atacama Cosmology Telescope (ACT). ACT's observing strategy and map making 
procedure allows an unbiased reconstruction of the modes in the maps over a wide range of multipoles. 
By directly matching the ACT maps to WMAP observations in the multipole range of 400 < I < 
1000, we determine the absolute calibration with an uncertainty of 2% in temperature. The precise 
measurement of the calibration error directly impacts the uncertainties in the cosmological parameters 
estimated from the ACT power spectra. We also present a combined map based on ACT and WMAP 
data that has high signal-to-noise over a wide range of multipoles. 
Subject headings: cosmology: cosmic microwave background, cosmology: observations 
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1. INTRODUCTION 



Observations of the cosmic microwave background 
(CMB) are one of the key probes of cosmology. 
Current balloon and ground based measurements 
(iSavers et al.1 120091: iBrown et alj l2 00l iFriedman et al l 



2001 iReichardt et al.1 l2009allbl; JSharp et al.1 120101:^ 
Dawson et al.l 20061: Sievers et al.1 120091 : iLueker et all 



.20101: IFowler et all I201(F and space-based experiments 
(jKomatsu et al.l 120101 ) measure the amplitude of cosmic 
microwave background fluctuations over a wide range 
of angular scales. All-sky maps of the Wilkinson 
Microwave Anisotropy Probe (WMAP) provide a large- 
scale measurement of the cosmic microwav e background 
with only a 0.2% calibration uncertainty (jJarosik et al.l 
120101 ). The scientific value of the small-scale ground 
and balloon-based measurements are maximized if these 
precise small-scale maps can be accurately calibrated to 
WMAP. 

Starti ng with the compa risons between FIRS and 
COBE (IGanga et al.l fl993l) . cross-calibration studies 
have both confirmed detections and provided com- 
mon calibration between experiments. General meth- 
ods for c ross-calibrati o n hav e been developed for CMB 
studies (|Ganga et al.l 119931 : iTegmarkl 
plied t o many different experiments: 



19991) and ap 



Reichardt et al.1 



(|2009bf ). (ACBAR with WM AP). IBrown et al.l (|2009l) . 
(QUaD with BOO MERang) iChiang et al.l <|2009T ). (BI- 
CEP with WM APVlXu et all (1200 lft. CS A SK and QMAP 
with CO BE'). INetterfield et al.l (I1997D. (SASK with 
MSAM). lAbroe et al.l (12004ft /Maxima with WMAP), 
iHernandez-Monteagudo eTall (12001 . (ARCHEOPS and 
WMAP). iMasi et al.1 " p006ft . (BOOMERanG and 
WMAP). An alternative approach to calibration relies 
on measuring known point sources, particularly planets. 
This ap proach has also b een used by many C MB exper- 
iments (jCrill et al.l 120031: iMason et all 12003ft and is ap- 
plied in section 12.21 

We present a new method based on cross-correlations 
with WMAP to accurately measure the absolute cali- 
bration uncertainty for ACT observations. The method 
we present here expands on earlier work in that we 
take full advantage of working in Fourier space. The 
Fourier space algorithms introduced in this paper are 
faster than the real-space algorithms. They are sim- 
ple and scale linearly with the number of pixels in 
the maps. Filtering noisy modes, deconvolving the 
beam and estimating the noise model is faster and eas- 
ier in Fourier space. The ACT data used in this pa- 
per were collected at 148 GHz and 218 GHz during the 
2008 observing season. The maps are solved for iter- 
atively using a preconditioned conjugate gradient code 
to obtain the unbiased maximum likelihood solution. 
For an overview of the ACT data re d uction pipeline, 
we re fer the reader to IFowler et al.l (|2010ft : iDunnerl 
(|2009ft . For the power spectrum results and cosmolog- 
ical parameters, see t he co mpanion papers (|Das et al.1 
|20T0[ iDunklev et all 12010ft . an d for clusters using 
the Suny aev-Zel'dovich e ffect dSunvaev fc Zel'dovich 



see iMarriage et al.l (|2010ft and iMenanteau et al 



For information on the t elesco p e facility, see 
IFowler et all (12007ft: iHincks et al.l (|2008l ): iSwitzer et all 
(|2008ft : iSwetz et al.l (20101 ). Details of the cryo- 
genic receiver and bolometric detectors are provided 




WMAP 



Fig. 1.— The ACT and WMAP maps used in this paper. The 
red region in the WMAP map shows the overlap between WMAP 
and the ACT data. 

inlNiemackl (12006ft: IMarriage et all (12006ft: IFowler et al.1 



(12007ft: iBattistelli et al.1 (120081): iNiemack et all (12008ft 
ISwetz et al.l (|2008f ): iThornton et all (|2008l ): iZhao et all 
(|2008ft . 

In Section [2] we briefly describe the data sets used for 
this analysis. In Section [3] we present a cross-correlation 
method to test the relative alignment of the ACT maps 
versus WMAP maps. Our calibration method is dis- 
cussed in Section [4] Section [5] presents our combined 
high resolution CMB map that is signal-to-noise domi- 
nated over a wide range of multipoles (I < 5000). 

2. DATA 

The cross-correlation method used in this paper re- 
quires a single map from WMAP and four maps with in- 
dependent noise properties from ACT at each frequency. 
Below we give a brief description of the three data sets 
that are used in this paper. Throughout this paper we 
use a to denote ACT data and w to denote WMAP data. 

2.1. ACT Data 

One of the strengths of ACT is its scanning strategy. 
The geographical location of the ACT telescope (the At- 
acama desert) enables cross-linking of observations. Ev- 
ery observing region is scanned along two different di- 
rections at constant elevation each night. Azimuth-only 
scans observe a sky field at the same elevation twice a 
night: once as the field rises; once as it sets. Sky rotation 
changes the scan angle, resulting in cross-linking, which 
allows (in principle) the reconstruction of all the modes 
in the map. The modes that are lost due to the corre- 
lated noise along the scan direction in the rising maps 
are won back using the setting maps and vice versa. 
Such cross-linking and the unbiased map-making algo- 
rithm used to make ACT maps yield a successful recon- 
struction of the modes in the maps down to multipoles 
of ~ 300. The combination of ACT's scan strategy and 
its unbiased map-making scheme enables a direct cali- 
bration to WMAP. However, since WMAP is noisy on 
small-scales and ACT's maps are dominated by atmo- 
spheric noise and poorly measured modes on its largest 
scales, our method compares the maps as a function of 
multipole over the 400 < I < 1000 range. 

The ACT maps used in this paper arc identical to the 
m aps used in IDas et alj (|2010ft . The area includes that 
m IFowler et al.l (|2010l ) and represents the highest sen- 
sitivity data from 2008 (148 GHz and 218 GHz). The 
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map resoluti on is 1.4' and 1.0' at 148 GHz and 218 GHz. 
respectively (|Hincks et al.l [20091 ) . The map projection 
used is cylindrical equal area with square pixels, 0.5' on 
a side. We divide the data into 12 patches. Each patch 
is 5° x 5° in size, and together they cover a rectangular 
area of the map from a = 00 h 22 m to 06 h 52 m (5?5 to 103°) 
in right ascension and from 5 = —55° to —50° in decli- 
nation. We divide our data-set into four equal subsets 
in time, such that the four independent maps generated 
from these subsets cover the same area and have sim- 
ilar depths. Therefore, for each patch described above 
we have four maps, each representing roughly a quarter 
of the time spent on that patch. We call these "season 
maps" . 

2.2. Calibration to Planets 

Observations of Uranus provide the initial calibration 
of the maps. Uranus was observed by ACT every few 
days during the 2008 season, yielding approximately 30 
usable observations. The time-ordered data from each 
observation are calibrated to detector power units and a 
map is produced. From each map, the peak response of 
the planet is determined, corrected for the temperature 
dilution due to the finite instrumental beam size (i.e. 
by multiplying by the ratio of Uranus's solid angle to 
the instrumental beam solid angle) and then compared 
to the Uranus temperature (converted to CMB differ- 
ential units) at the band effective frequency. The cali- 
bration factors from detector power units to CMB dif- 
ferential temperature are compared to precipitable wa- 
ter vapor (PWV) measurements to fit a model of atmo- 
spheric opacity, and the season map calibration factor is 
obtained by evaluating the fit at the season mean PWV 
of 0.5. The uncertainty in the Uranus-based calibration 
is 7%, and is dominated by the 5% uncertainty in the 
planet's temperature. 

We take the brightness tcmperaturifl of Uranus to be 
107 ± 6 K and 96 ± 6 K for the 148 GHz and 218 GHz 
bands respectively. These temperature s are based on a 
reproc essing of the data presented by iGriffin &: Ortonl 
(|1993|) . in combination with WMAP 7- year measure- 
ment s of Mars and Uranus brightnesses (Weiland et al. 
l20T0t) . Griffin and Orton measured flux ratios between 
Uranus and Mars at several wavelengths from 350 Lim to 
1.9 mm. These were calibrated to abs olute units u sing 
an extrapolation of the Wright model ()Wrightlll976[ ) for 
Mars temperature, and the resulting Uranus brightness 
temperatures were fit with a third-order polynomial in 
the logarithm of the wavelength. We have followed the 
same procedure, but with two modifications. The first is 
that the WMAP measurements of Uranus temperature 
at 94, 61 and 41 GHz have been included in the fit. The 
second is that the WMAP comparison of Mars temper- 
ature at 94 GHz (3.2 mm) to the Wright model at 350 
|xm is used to pin the long wavelength end of the Mars 
brightness extrapolation. The resulting Uranus bright- 
ness temperatures in the two bands are 5% and 3% lower 
than those obtained by Griffin and Orton. This differ- 
ence is due almost entirely to the recalibration of the 

1 Throughout, the brightness temperature refers to the temper- 
ature of a black-body having the observed in-band spectral radi- 
ance. The equivalent Rayleigh- Jeans temperatures of Uranus for 
our bands arc 103.5 and 90.7 K. 



TABLE 1 

THE RMS NOISE FOR EACH WMAP DIFFERENCING ASSEMBLY. 





Wi 


W 2 


w 3 




(Hv(niK) 


5.906 


6.572 


6.941 


6.778 



Mars model. 



2.3. WMAP Data 



The WMAP map uses co-added inverse-noise weighted 
data from seven single year maps and four differenc- 
ing assemblies at 94 GHz0 The maps are foreground 
cl eaned (using th e foreground template model discussed 
in iHinshaw et al.l <|2007f )) and are at HE ALPix resolu- 
tion 10 (iVsido = 1024), with 3.5' pixels. Single year maps 
are multiplied by the pixel weights based on pixel noise 
evaluated with the expression 



a 2 (n) 



Mr 



/N obs , 



(1) 



where N b s is the number of observations at each map 
pixel which is directly proportional to the statistical 
weight and aw is the noise for each differencing assem- 
bly. The four noise factors in mK units are given in Table 

m 

We cut a rectangular region from the resulting 
HEALPix map corresponding to the ACT southern strip 
and project it onto the CEA coordinates. Both the 
WMAP and ACT data are cut into the same 12 patches. 
Figure [2] shows one of these patches and the correspond- 
ing patch from the ACT measurements. While both 
maps shows similar hot and cold spots, the WMAP beam 
has smoothed out small scale structure and the ACT 
maps contain large-scale atmospheric noise. 

Because the ACT maps have poorly measured modes 
on the largest scales, the WMAP and ACT maps are 
filtered by a high-pass filter F c (l) in Fourier space be- 
fore being cut into 12 patches. The high-pass filter is a 
smooth sine-squared function in Fourier space given by 

F c (l) = sin 2 x(0©a-/ mi „)©(Z m ax-0+ea-'max), (2) 

where x(l) = (tt/2)(Z - J mm )/Gmax - 'mm) and O is the 
Heavyside function. We choose l m - m = 100 and Z max = 
500. The final power spectrum is corrected for this filter 
as well as for the effects of the beam and pixel window 
functions. 

3. COMPARING MAPS 

We test the alignment of the maps by computing their 
real-space cross-correlation function as a function of off- 
set x. This is done by inverse Fourier transforming the 
2-dimcnsional cross-correlation in Fourier space: 

M aXw (x) ='^2a(l)w*(l)exp(il-x). (3) 
I 

Here a(l) and w(l) are the Fourier transforms of the 
ACT and WMAP maps respectively, and the tilde de- 
notes high-pass filtering of those maps, a(l) = a(l)F c (l) 
and w(l) = w(l)F c (l). The high-pass filtering is impor- 
tant here because otherwise long wavelength noisy modes 

2 from http : //lambda . gsf c . nasa . gov 
,? http: //healpix . jpl .nasa . gov | 
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Fig. 2. — One of the 12 patches used for cross-correlations from ACT 148 GHz (right), and WMAP data (left). Long wavelength modes 
have been removed from both maps using the high- pass filter described in eqn. J2). The ACT map has been convolved with the WMAP 
W-band beam. The two maps represent the same area in the sky, observed by two experiments. The similarity between structures in 
the two maps is the key feature that we use for our cross-correlation studies. The ACT maps have large-scale atmospheric noise and the 
WMAP maps are dominated by detector noise on the smallest scales. These two types of noise are completely uncorrelated. 



dominate the cross correlation and hide the true signal. 
The resulting 2D cross-correlation function, M a 

x w ) can 

be written as 

M axw (x) = C(x) ® B{x) + n{x), (4) 

where B(x) is an effective beam between the two maps, 
C{x) is the correlation function of the two sky maps and 
n(x) is the noise fluctuations in the cross-correlation that 
gets smaller as we add more area. For a perfectly aligned 
pair of maps, the 2D cross-correlation function, M axw (x) 
peaks at x = and its width is determined by the correla- 
tion length of the field and the beams. For our maps, the 
result of this is shown in Figure [3J The cross-correlation 
function is cleanly peaked at zero lag, which shows the 
maps are aligned and correlat ed. The noiseless sky sim- 
ulations of lSehgal et al.l (|2010f ) provide a confirmation of 
these results. We convolve the simulated maps with the 
WMAP and ACT beams to generate simulated noiseless 
maps of WMAP and ACT data respectively. The maps 
are then filtered with the high-pass filter of eqn. © and 
cut into 12 patches. The real-space cross-correlation of 
the simulated maps is computed using eqn. (j4]). The 
result is plotted in Figure [3] (dashed blue lines) . The 
agreement between cross-correlations computed from the 
data and the simulations is striking. We fit a Gaussian 
curve to the radial profiles of the 2D cross-correlations. 
At small scales (< 10') a a = 7.7' Gaussian profile is 
a good fit, and the result of the simulations agrees well 
with the data. 

4. CALIBRATION USING CROSS CORRELATIONS 

This paper uses four kinds of cross-power spec- 
tra, ACT x WMAP, ACT x ACT ACTi 48 xACT 2 i 8 , and 
WMAP x WMAP to provide CMB-bascd calibrations of 



the ACT maps. The advantage of working with the cross- 
power spectrum is that the noise in the two maps is in- 
dependent and therefore uncorrelated. Hence a precise 
modeling of the noise is not needed. 

If the noise in the two maps is uncorrelated, the cross- 
power spectrum provides an unbiased estimator of the 
underlying power spectrum. Cross-correlations can be 
done between pairs of maps from two different experi- 
ments. Because the two experiments have different sys- 
tcmatics that do not correlate, the cross-spectrum pro- 
vides a good estimate of the true underlying power spec- 
trum that is nearly independent of the noise properties 
in the maps. 

4.1. Power Spectrum Method 

The power spectrum method we use is the Ad aptive 
Multi-Taper Method (AMTM) of lDas et al.l (|200l . Also 
we would like to have the maximum resolution possible 
in Fourier space to have as many bins in the power spec- 
trum as possible, in order to maximize the number of 
independent measurements of the calibration factor as 
it is described below. The size of the bins is limited by 
the fundamental frequency determined by the smallest 
side of the maps. For the maps we are using, this is 
SI = 72. Using the AMTM with multiple tapers would 
result in smaller errors at the large I regime of the power 
spectrum, but it increases the size of the independent 
bins. For this reason we use AMTM with one taper at 
resolutio n Nres = 1 and thus no iterations are neces- 
sary (see iDas et al.l ()2009l ) for the details of the AMTM 
method). We use slightly larger bin sizes than the fun- 
damental frequency (SI = 90) to guarantee that the bins 
are uncorrelated. This is further tested and verified by 
Monte Carlo simulations. The prewhitening method of 
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Fig. 3. — Comparing ACT and WMAP data sets: Two-dimensional cross-correlations of WMAP and ACT maps in real-space (inset 
panels) and their one dimensional radial profiles (solid red curves) for the 148 GHz maps (left) and 218 GHz maps (right). 2D cross- 
correlations peak at zero lag. This solidifie s that the maps are aligned and the pointing agrees on average over the whole map area. Results 
obtained from the noiseless simulations of lSehgal et al. (2010) are also plotted (dashed blue lines). The data correlation functions agree 
with the simulations. A a = 7.7' Gaussian profile is a good fit to the 2D correlations at small scales (< 10'). 



iDas et ail (|2009() is designed to reduce the dynamic range 
of the Fourier components of the maps. This is impor- 
tant when one is interested in measuring the damping 
tail of the CMB (I > 1000) where the slope of the power 
spectrum is steep. However we do not use prewhiten- 
ing as we are working in the mildly colored regime of 
the power spectrum, I < 1000 where the dynamic range 
of the power spectrum is not large. Using prcwhitcning 
does not affect the power spectrum at the I ranges of our 
interest. 

The ACT power spectrum (both for 148 GHz and 
218 GHz maps) is computed using cross correlations of 
the four season maps. The power spectrum of each patch 
is the average of the six cross-power spectra: 



patch: 



C\ = 



12 



l</3;<4 

E 

:;,ft;a,<ft 



c\ 



aiPi 



(7) 



The ACT x WMAP cross-power spectrum, Cf w , is 
given by the average of the 12 -patch cross-power spec- 
tra. The spectrum in each patch is given by the average 
of the four cross-spectra between the WMAP map and 
each of the four ACT season maps for that patch, 



c\ 



1 



E c r 



(8) 



Ci 



l</3i<4 

E 

i,Pi;cii<Pi 



(5) 



where a, indexes the four season maps of the ith patch 
in ACT data and 



d6a ai (l)w*(l), 



(9) 



where C" 1 ^ 1 arc cross-power spectra of season maps for 
the patch i, and a» and pi index the four season maps of 
that patch. The final ACT power spectrum is given by 
the average of the 12 patches 



Ci 



1 

12 



12 

E 

i=i 



ci 



(6) 



Each of these 12 power spectra is an independent mea- 
surement of the ACT power spectrum. We use the vari- 
ance from the 12 power spectra values at each I bin as 
a measure of the error on the power spectrum. This 
metho d agrees well with the analytical estimate of the 
errors (jFowler et al.ll2010D . 

The ACTi4g x ACT218 cross-power spectrum is com- 
puted in a similar way, but in this case cti and /3i cor- 
respond to the season maps from 148 GHz and 218 GHz 
data respectively, and there are 12 cross-spectra for each 



where Ni is the number of modes in each bin in Fourier 
space. 

The window function deconvolution of the power spec- 
tra is done on the average cross-power spectrum of each 
patch. We have tested our pipeline on 1000 Monte Carlo 
simulations to confirm that the power spectrum method 
we use (including the 2D noise weighting and the /-space 
masking) is unbiased and that the covariance matrix of 
the po wer spectrum bin s is diagonal. The patch sizes 
used in IDas et ail (|2010l ) are three times larger than the 
pa tch sizes u s ed in this paper, but the number of patches 
in IDas et al.l ([2010) is three times fewer than that used 
here. So the total areas covered by both cutting methods 
are identical and for the same binning, the final power 
spectra obtained from the two methods agree with each 
other to better than 1% fractional error. 

The straight binning of the two-dimensional power 
spectra is the simplest but not the best. Down- weighting 
noisy regions in the 2D power spectrum space before bin- 
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Fig. 4. — Average of 72 noise estimates (\xK ) in 2-dimensions for 148 GHz (left) and 218 GHz (right). Each estimate is computed from 
the auto-spectrum of a difference map obtained from differencing two quarter-season maps from the same region of the sky. The sky signal 
cancels in the difference maps and only noise remains. The anisotropy of the noise power spectrum is the reason we use the se noise models 
to dow n-weight the noisy regions in 2-dimensional power spectra before binning. These weights arc different from those of Marriage ct al. 
[|2l)T0h in that the CMB is the main signal for us and it docs not contribute to the noise model. 



ning is a useful technique that is adopted to improve the 
power spectrum method. Below we discuss our method 
for estimating the noise model and the 2D noise weight- 
ing along with the /-space masking. 

4.2. Noise Model 

The season maps for each patch share the same sig- 
nal but have independent noise properties. Differenc- 
ing the two season maps removes the sky signal and 
leaves behind a linear combination of the noise in the 
two maps. The resulting map, which we call the "differ- 
ence map," is a noise map and its auto-power spectrum 
can be used for estimating the noise model for that patch 
(see lMarriage et al.l (|2010l ) for a more detailed discussion 
on noise weighting in Fourier space) . Figure [4] shows the 
stacked noise models given by the average of 72 two di- 
mensional power spectra obtained from 6 difference maps 
per patch, for 12 patches. 

4.3. Noise Weighting 

The ACT and WMAP maps have very different noise 
properties. The WMAP noise spectrum is dominated 
by the detector noise and is nearly constant over a wide 
range of the wave- vector /, whereas ACT noise has more 
detailed structure. The ACT noise is dominated by at- 
mospheric noise on large angular scales, and by the de- 
tector noise on small scales. Some directions in the two- 
dimensional ACT power spectrum are more noisy than 
others. We use our best estimate of the noise for each 
map to down weight the noisy parts of the spectra before 
angle-averaging the two-dimensional spectra. The noise 
for the ACT- ACT cross-spectrum is N% a (l) and for the 
ACT- WMAP cross-spectrum it is N aa (l)N ww (l), where 
N aa (l) is the ACT noise spectrum in two dimensions 
as shown in Figure 0] and N ww (l) is the WMAP noise. 
Noise weighted power spectra are then angle-averaged 
and binned in /. In the end, the spectra are divided by 
the beam and relevant pixel window functions. If we 
denote the two dimensional spectra by P(l), the final 
ACT- ACT cross-spectrum is obtained by 



C aa (l) 

where b a (l) is the 

(jHincks et al.l [2009h . 



fdOP aa (l)/N2 aa (l) 
F?(l)bl(l)f<10/K a (iy 



(10) 



ACT beam in Fourier space 
F c (l) is the high-pass filter in 



Fourier space and J d0 represents angle-averaging and 
binning in /. We do not correct for the ACT pixel win- 
dow function as it is close to unity in the range of interest, 
I < 1000. 

The ACT- WMAP cross-spectrum is obtained from 



C aw (l) 



J d9P aw (l)/N aa (l) 



F*(l)b a (l)b w (l) Pw (l) f d6/N aa (iy 



(11) 



where b w (l) is t he W band beam of WMAP 
(jJarosik et al.|[2010h and p w {l) is the pixel window func- 
tion of the resolution 10 HEALPix maps corresponding 
to N side = 1024L with 3.5' pixels. We use HEALPix in 
Python (healpjp) to compute the pixel window function 
at the integer harmonic indices and then bin it in I to 
get p w (l)- The WMAP noise term gets canceled in the 
above noise weighting as it is constant. 

We mask a vertical band of width Al x = 180 in the 
power spectra before angle-averaging them. This makes 
sure that our estimates of the power spectra are not af- 
fected by the striping effects that are present at a narrow 
band around / = 0. For this analysis, we use a /-space 
mask at l x = [— 90, 90], For a detailed d escription of the 
/-space masking see iFowler et all (|20Tfjt) . 

4.4. Calibration: ACT 148 GHz 

We assume that the calibration factor is constant on 
various scales and that the ACT, a(x), and WMAP , 
w(x), maps can be represented as 

a(x)=a- 1 AT sky {x)®B a (x)+N a (x), (12) 
w(x) = AT sky (x) ® B w (x) + N w (x), 

where a is the calibration factor ([Dunnerl 120091) . 
AT s k y (x) is the sky temperature signal, N a and N w arc 
ACT and WMAP noises respectively, B a is the ACT 
beam and B w is the WMAP beam. The calibration fac- 
tor a appears in the cross and auto-spectra of ACT, and 
it can be estimated through relevant ratios of the auto- 
and cross-spectra. We define 



ai (i)=cr/cr 

a 2 (l) = C l /Gj 



(13) 



' http : //code . google . com/p/healpy 
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in which Cj uw is the average 94 GHz (W-band) power 
spectrum of the WMAP 7-year data (jLarson et al.ll201(# l 
binned in the same way as the ACT auto- and cross- 
spectra. The above three measures of the calibration 
(cqn. (|13p ) are not independent, but they can be used 
together to test the consistency of our results over a wide 
range of angular scales. Variation of a versus I is a sign of 
scale dependent calibration factor. For the ACT 148 GHz 
maps, the calibration factor is constant in the range of 
400 < I < 1000, and different measures are consistent 
with each other within the errors. We restrict our analy- 
sis to this region. For I < 300 the ACT data is dominated 
by the atmospheric noise. Beyond I > 1000 WMAP maps 
arc resolution limited. 

The power spectrum bins are chosen such that the co- 
variance between bins is negligible. Therefore the ra- 
tios of the spectra, a(l), provide independent measure- 
ments of the calibration factor at every I bin. For the 
ACT 148 GHz maps the a(l) are flat over the range of 
400 < I < 1000. Hence the overall calibration factor a 
can be estimated from the average of the a(l) values in 
that I range. Since the covariance between the a(l) values 
at different I is negligible, the variance of the quantities 
that are used in averaging is a good measure of the error 
on the mean. Therefore we obtain: 

Ql = 1.01 ±1.9%, (14) 
a 2 = 1.00 ±2.1%, 
d< 3 = 1.00 ± 1.4%. 

The uncertainties quoted are the error on the mean 
derived from the variance of the a(l) for every case. We 
use the Anderson-Darling statistic to test the normality 
of the a(l) measures. The A 2 statistic that quantifies de- 
viations from normality in this test becomes large when 
data points deviate from normality. If A 2 > 0.751, the 
hypothesis of normality is rejected at the 95% confidence 
level (for a 90% confidence level it is 0.632). We compute 
A 2 for the three measures of the calibration factor defined 
above. The result is A 2 ( ai ) = 0.27, A 2 (a 2 ) = 0.26 and 
A 2 (a 3 ) = 0.43. Therefore at the >10% level, the data 
used in computing the average calibration factor do come 
from a normal distribution. 

Among the above three measures of the calibration fac- 
tor, &3 has the smallest error. The reason is that it uses 
the sky averaged WMAP power spectrum as the estimate 
of the WMAP spectrum. Smaller errors in this quantity 
translate into the smaller error in the a^. However all 
three measurements are in agreement with each other. 
For comparison, the Uranus calibration is 0.99 ± 7.0%. 
We use 2% as the calibration error for the 148 GHz maps. 

The calibrated power spectra are shown in Figure 
The ACT 148 GHz - WMAP cross-spectrum is shown 
with blue filled circles, and single frequency spectrum of 
the ACT 148 GHz maps is shown with red triangles. The 
average 94 GHz (W-band) power spectrum of the WMAP 
data binned in the same way as other spectra are binned 
(light green boxes) and the theory power spectrum based 
on WMAP7 best fit parameters (solid black line) are also 

5 http: //lambda. gsfc .nasa. gov/product/map/dr3 



plotted for reference. ACT x ACT power spectrum has 
large uncertainties on large scales due to the large scale 
noise in the ACT maps. WMAP power-spectrum has 
larger errors on smaller scales where WMAP detector 
noise dominates and maps are resolution limited. 

4.5. Calibration: ACT 218 GHz 

We measure the calibration factor for the 218 GHz 
maps in a similar way. The result is a calibration factor 
with ~ 7% fractional uncertainty in temperature. The 
larger uncertainty on the ACT 218 GHz calibration fac- 
tor is due to the higher noise level in those maps than the 
148 GHz maps on large angular scales. The calibration 
factors determined from cross-correlations with WMAP 
and from Uranus observations are consistent to ~ 7%. 
We use 7% calibration error, based on Uranus, for the 
218 GHz maps. 

4.6. Calibration systematics 
Several factors can systematically affect our results: 

• Stability vs. time and space: In order to examine 
the stability of our calibration results in time and 
space, we compare the average calibration factors 
in the 12 patches. The average ai in every patch 
is computed by averaging the calibration factor for 
different I values in that patch. The scatter of the 
average calibration factors in 12 patches gives us a 
measure of stability of our results versus the posi- 
tion in the sky and the time that the patch was ob- 
served. We verify that the patch-to-patch variation 
of the calibration factor is consistent with the la 
spread in the average calibration factor. We have 
also checked the independence of the calibration 
factor on the area of the map. The calibration fac- 
tor m easured on the 228 deg 2 data of iFowler et al.l 
( 2010 ) is the same as that measured on the ^300 
deg map used in this work. 

• Beam uncertainty and w indow func t ion n ormal- 
ization: As discussed in iPage et aD ()2003l ). lim- 
ited knowledge of the beam profiles leads to un- 
certainties in the experimental window function. 
These uncertainties result in a distortion of the 
power spectrum, which in turn affects our estima- 
tion of the calibration factor. One of the uncertain- 
ties is the uncertainty in the normalization of the 
ACT window function, which appears as an overall 
change in the calibration factor. The beam transfer 
function we use is normalized to unity at I = 700 
and thus has no uncertainty at that The uncer- 
tainty in the beam transfer function is taken into 
account when estima ting cosmological parameters 
(jDunklev et all 12010( 1 and we do not deal with it 
separately here. 

• Pointing reconstruction error: The err or due to the 
pointi ng reconstruction is discussed in iJones et al.l 
(|2006[ ) and causes a correlated distortion of the 
power spectrum at different bins. Absolute de- 
tector array pointings for the ACT data are es- 
tablished with 3.5" precision through an iterative 
process in which the absolute pointing is adjusted 
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Fig. 5. — ACTxWMAP cross-power spectra. The ACT 148 GHz - WMAP cross-spectrum is shown with blue dots and the single 
frequency spectrum of the ACT 148 GHz maps is shown with red triangles. The average 94 GHz (W-band) power spectrum of the WMAP 
data (light green squares) and the theory power spectrum based on WMAP7 best fit parameters (solid black line) are also plotted for 
reference. The gray band shows the cosmic variance. 



based on offsets of ACT-observed radio source lo- 
cations with respect to source locations taken from 
the Australia Telesc ope 20 GHz (AT20G) survey 
(jMurphv et al.ll2009D . The bias induced by this er- 
ror in the window function is not significant at the 
scales of our interest. 

• The band center uncertainty: Our calibration 
method is based on the CMB cross correlations in 
two different frequencies. The CMB spectrum in 
the frequency range of this work is flat and there- 
fore uncertainty in the frequency band center of the 
experiment will not affect our results. 

5. COMBINING MAPS 

The ACT and WMAP views of the CMB sky are 
complementary. WMAP accurately measures large-scale 
CMB features and ACT provides a high resolution im- 
age of small scale features. As we have shown in sec- 
tions 3 and 4, these maps are consistent and now cross- 
calibrated so that they can be combined to make a CMB 
map that is signal-to-noise dominated over a wide range 
of scales. The WM AP data are signa l dominated on 
large scales, I < 548 (|Larson et al.ll2010l) and the detec- 
tor noise dominates at smaller scales. The ACT data are 
dominated by the unconstrained modes on large scales 
that arise from a combination of instrument properties, 
scan strategy, and atmospheric contamination. The ACT 
data have a good signal-to-noise ratio on intermediate 
scales, 500 < I < 5000, and become noise dominated 
beyond that. The basics of our method for combining 
WMAP a nd ACT is simila r to the real-space method de- 
scribed in iTegmarkl (|1999f ) . Working in real space has 
the disadvantage of having to deal with large matrices 
for large maps. Instead, we work in Fourier space. We 
construct a linear combination of the two maps in Fourier 
space by inverse noise weighting them such that the less 
noisy map dominates at each scale (these maps are not 



high-pass filtered). The WMAP data has lower angu- 
lar resolution than ACT and the pixel size of WMAP is 
much bigger than that of ACT as well. We first decon- 
volve the WMAP data to the ACT angular resolution in 
Fourier space. Then we use inverse variance weighting 
to combine the two maps. The construction is given by 



M(l) 



,(l)a(l) 



erf, 



,{iy 



(15) 



where a aa and a ww are the ACT and WMAP noise spec- 
tra in two dimensional Fourier space. The above com- 
bination picks up WMAP contributions on large scales 
and ACT contributions on small scales and is signal dom- 
inated over a large range of I. The combined Fourier 
map, M(l), is transformed back to real space. The ACT 
noise model, a aa , is computed in the same way as de- 
scribed above. WMAP noise is estimated by a white- 
noise, i.e. constant, before beam deconvolution. For 
beam deconvolved WMAP maps, the noise model, cr ww , 
is the WMAP noise multiplied by the two dimensional 
beam window function in Fourier space. The resulting 
map is shown in Figure [6l This map encompasses all 
components of a high resolution CMB map (CMB on 
large scales, point sources and clusters on small scales) 
in one single map. 

6. CONCLUSION 

We have presented a method for calibrating ACT maps 
using cross-correlations with WMAP. The fractional un- 
certainty in the temperature calibration factor obtained 
in this paper is ~ 2% and the absolute calibration agrees 
well with that derived from Uranus observations. The 
calibration uncertainty is inversely proportional to the 
total area of the maps used for cross-correlations. Adding 
more area will decrease the uncertainty in the calibration 
factor. We have also presented fast methods for compar- 
ing and combining CMB data sets. Using these methods, 
we have combined ACT with WMAP data to construct a 
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Fig. 6. — A 84 deg 2 combined map obtained by combining WMAP and ACT 148 GHz maps. This high resolution map encompasses all 
components of the microwave sky: the large scale structure is the CMB, hot (red) dots are point sources and some of the small cold (blue) 
spots are clusters of galaxies. 



high resolution map that has a good signal-to-noise ratio 
over a wide range of angular scales. 
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